Electronic transport and quantum localization effects in organic semiconductors. 
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We explore the charge transport mechanism in organic semiconductors based on a model that 
accounts for the thermal intermolecular disorder at work in pure crystalline compounds, as well 
as extrinsic sources of disorder that are present in current experimental devices. Starting from 
the Kubo formula, we describe a theoretical framework that relates the time-dependent quantum 
dynamics of electrons to the frequency-dependent conductivity. The electron mobility is then calcu- 
lated through a relaxation time approximation that accounts for quantum localization corrections 
beyond Boltzmann theory, and allows us to efficiently address the interplay between highly con- 
ducting states in the band range and localized states induced by disorder in the band tails. The 
emergence of a "transient localization" phenomenon is shown to be a general feature of organic 
semiconductors, which is compatible with the bandlike temperature dependence of the mobility ob- 
served in pure compounds. Carrier trapping by extrinsic disorder causes a crossover to a thermally 
activated behavior at low temperature, which is progressively suppressed upon increasing the carrier 
concentration, as is commonly observed in organic field-effect transistors. Our results establish a 
direct connection between the localization of the electronic states and their conductive properties, 
formalizing phenomenological considerations that are commonly used in the literature. 



I. INTRODUCTION 



Remarkable progress has been made in recent years 
in understanding and improving electronic transport in 
organic semiconductors (OSC) and devices. Mobilities 
exceeding IQcm? /V s are now measured in an increasing 
number of organic field-effect transistors (OFETs) based 
on single crystalJ^^. Such values are orders of magnitude 
lower than those attainable in inorganic semiconductors, 
and are indicative of extremely short electronic mean-f ree 
paths — on the order of the inter-molecular distance^^ 
— causing a breakdown of the basic assumptions under- 
lying band transport. This occurs despite the relatively 
modest coupling of the carriers with intra-molecular vi- 
brations, whi ch rules out the presence of polarons in such 
materials.!^ It is currently believed that the mobility 
in crystalline organic semiconductors is intrinsically lim- 
ited by the presence of large thermal molecular motions, 
which are a direct consequence of the weak Van der Waals 
inter-molecular bonds .'^^'S! Deviations from the perfect 
crystalline arrangement act as a dynamical source of dis- 
order on the already narrow electronic bands arising from 
the TT-intermolecular overlaps, inducing a localization of 
the electronic wavefunctions on the timescale of the inter- 
molecular vibrations themselve^^ — a phenomenon that 
is not described by the semiclassical Boltzmann theory of 
electron-phonon scattering, nor by the classical Marcus 
electron transfer theory. Theorie s based on such "tran- 
sient electron localization"'SniIlll are able to explain the 
power-law decrease of the mobility with temperature ob- 



served in ultrapure organic semiconductors^^ as well as 
the optical conductivity data available in OFET ^ I ^^ I ^^ l 

Experimentally, the intrinsic mobility of organic semi- 
conductors is still difficult to observe in practical OFET 
devices. Even when polaronic self-trappin^i^ and dipo- 
lar disordei^ induced by the interface polarizability are 
avoided by using non-polar gate dielectrics or suspended 
samples, the carrier mobility in OFETs is still affected 
by extrinsic sources of disorder, related to the presence 
of structural defects or to the interface roughness. Ex- 
trinsic disorder favors the formation of trapped states in 
the band tails , at energies located below those of the in- 
trinsic carrierPESl. As a result, depending on the device 
quality, a crossover from an intrinsic transport regime to 
a thermally activated (tr appe d) regime is observed upon 
lowering the temper at ure^^El^ or the intrinsic regime can 
be completely washed out if the disorder is sufficiently 
strong as occurs in polycrystalline films.'^il 

As is clear from the above discussion, a proper descrip- 
tion of the transport mechanism in both pure and disor- 
dered OSC requires a method that (i) goes beyond both 
Boltzmann and Marcus approaches and (ii) is able to 
describe the interplay between highly conducting states 
in the band range and weakly mobile states induced by 
disorder. This is achieved here by applying a recently 
developed theory of charge transport based on the Kubo 
formula for the electrical conductivity combined with a 
suitable relaxation time approximation (RTA) on the 
current-current correlation function, which takes quan- 
tum localization effects into account .E^^^H xhe present 
method has already been successfully applied to analyze 
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the quantum transport properties of quasicrystalfl^Ull 
and to address the role of defects in graphene^^]. It has 
also been shown to provide an efficient description of the 
transient localization phenomenon in pure OSC!^, as it 
gives access to the time-resolved diffusivity and localiza- 
tion length of electronic states. By addressing the same 
quantities resolved in energy, we show here that this the- 
oretical framework also establishes a direct relationship 
between the existence of competing electronic states at 
different energy scales and the resulting transport prop- 
erties. Accordingly, both the intrinsic transport mecha- 
nism of clean organic semiconductors and the crossover to 
a thermally activated motion in the presence of extrinsic 
disorder are rationalized in terms of the relative weight 
played by strongly localized tail states and more mobile 
electronic states in the band range. The increase of the 
mobility observed in OFETs upon injecting a sufhciently 
large density of carriers is also naturally explained within 
this scenario. 

The paper is organized as follows. In Sec. |ll] we in- 
troduce the formalism relating the quantum diffusion of 
electrons to the Kubo response theory. Based on this 
formalism, in Sec. |III| we briefly describe the semiclassi- 
cal approximation used in Ref. |I2l and then derive the 
relaxation time approximation to be used here. A model 
relevant to organic semiconductors and devices is intro- 
duced in Sec. IIVI The results obtained in the limit of 
low carrier concentration are presented in Secti on [V] and 
their density-dependence is analyzed in Sec. |VI| The 
main conclusions are drawn in Sec. IVIIl 



II. GENERAL FORMALISM 

A formalism that relates the quantum diffusion of elec- 
trons, i.e. the quantum mechanical spread of the electron 
position with time, to the optical conductivity was intro- 
duced in Refs. l22l23l for metals and generalized to semi- 
conductors in Ref. [HI The main steps of the derivation 
are reviewed here. Readers not interested in formal de- 
velopments may skip this Section and move on directly 
to Section III. 



A. Optical conductivity and time-resolved 
diffusivity 

We start from the Kubo formula that relates the re- 
sponse of electrons to an oscillating electric field to the 
current-current correlation function (say, along a;)^^: 



1 



Re 



dte'("+^*>([J,(t),J,(0)]). (1) 



Here 5 is a small positive number enforcing convergence, 
n is the system volume, and we have set h = 1. Denoting 
the retarded current-current correlation function as 



(2) 



and its its Fourier transform as C_ 
Eq. ([T]) can be expressed as 



the Kubo formula 



Re C-{uj). 



(3) 



A relation between the mean square particle displace- 
ment and the current correlations can now be obtained 
through the retarded current-current anti-commutator 
correlation function, 



C+{t) = 9{t){{X{t),U0)}). 



(4) 



Writing the current operator in terms of the velocity 
operator, J — eV — edX/dt, and performing the time 
derivative we see that this function is directly related to 
the mean square displacement AX^{t) — (|A(t)— A(O)p) 
of the total position operator X{t) = J^i^ii^) along the 
chosen direction, (with the position operator for the 
i-th particle) via 



dAX^{t) 



dt 



C+{t')dt'. 



(5) 



This defines the instantaneous diffusivity of a system of 
N quantum particles. 



IdAX^it) 
2 dt ■ 



(6) 



Introducing the mean square displacement reached by 
the A^-particle system over a typical timescale r as 



dte'^' 



^dAX^t) 



dt 



(7) 



and using the properties of Laplace trasforms of deriva- 
tives, Eq. ([5]) yields the following relation between the 
mean square displacement and the Laplace transform 
C+(p) of the anti-commutator correlation function,'2Sl 



C+ip=l/T) = e- 



(8) 



The above equation shows that the quantity (7+ (p) has a 
precise physical meaning: it corresponds to the diffusivity 
of the electronic system averaged on a timescale r = p~^. 

Because the functions C+ and C_ are related by the 
detailed balance condition, which in Fourier space reads 



i?eC_ (uj) = tanh 



f3uj 



ReC+{uj) 



(9) 



(with /? the inverse temperature, see appendix [A| , the 
two relations Eqs. ^ and ([8| are are not independent. 
Indeed, by expressing the right-hand side of E q. pi in 
terms of the Laplace transform C+(p), via Eq. ( |A3[ ), we 
obtain an expression relating the mean square displace- 
ment L'^{t) and the optical conductivity cr(w): 



duj 



2uj 



t: + (1/t)2 tanh(/3w/2)' 



(10) 
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Remarkably, this relationship allows to address the time- 
resolved diffusion of electrons from the knowledge of the 
optical conductivity, which is a spectral property. An 
analogous equation was derived in Ref. for the in- 
stantaneous spread AX^(f). 



B. d.c. conductivity and mobility 

From the equivalence of the two formulations Eqs. ([3| 
and ([s]) we can derive a generalized Einstein relation 
connecting the electrical conductivity to the extensive 
diffusion coefhcient V, which is valid for quantum A'^- 
particle systems. By definition, a system is diffusive if 
the diffusivity at long times tends to a constant value, 
limt_j.oo T^it) = 2?. In the limit r — >■ cx) the integral in 
Eq. ([?]) is then dominated by such asymptotic diffusive 
behavior leading to 



V 



lim 

r— f oo 



LHt) C+ip = 0) 



2t 



2e2 



(11) 



Conversely, the above equation shows that reaching a 
finite localization length i(T oo) in the long time limit 
implies a vanishing diffusion coefficient. 

Using Eq. Q together with the definition of the d.c. 
conductivity from Eq. ([s]) as the limit 



lim 



(12) 



and observing that limp_j.o C+{p) = limi^^o C+{u!)/2 (see 
appendix \M we can write 



-V = 



lim 



(13) 



Our definition of the extensive diffusion coefficient T) 
for the N-particle system differs from the usual single 
particle diffusivity, which we denote D. The latter is 
an intensive quantity, defined as the ratio between the 
conductivity and the charge-charge susceptibility'^, so 
that 



with n the density and /i the chemical potential. The 
density-dependent proportionality factor on the r.h.s. 
corresponds to the number of particles that can actu- 
ally move, i.e. the compressibility times the (thermal) 
energy interval. Accordingly, the mobility of electrons 
can be defined at any finite density via 



eD 
k^ 



eksT 



d7i 



(15) 



C. Energy resolved quantities 



regions of the electronic spectrum: tail states generated 
by disorder below the band edges behave differently from 
states within the electronic band. It is therefore useful 
to decompose the response of the electronic system into 
contributions from states at different energy scales. '^^'^Il 
This can be done by exploiting the following expression 
of C_(a;) as an energy integral (see appendix [A|) : 



C_(w) = 7r J diy[f{iy) - f{uj + iy)]tr[p{iy)Jp{iy + oj)J], 

(16) 

where and f{v) = [e^'~''-f') + l]-i is the For mi function, 

= --Im^-^ (17) 

is the spectral operator from which the DOS p{i') can 
be obtained as p{i') — trp(iy), and H is the Hamiltonian 
operator. Defining 



B{iy)^tr[p{iy)Jp{i^)J], 



(18) 



the d.c. conductivity is readily obtained from Eq. ( 12 ) 

as 



dv 



(19) 



[note that there was a misprint in the definition of B{i') 
in Ref. J3J. We see from the above equation that the 
total conductivity of an electronic system arises from an 
average of B{v) over all electronic states, weighted by the 
corresponding statistical population. For example, in a 
system at finite electron density and low temperature, 
because the derivative of the Fermi function is peaked at 

~ /i, the conductivity is determined by the electrons in 
proximity (within ksT) of the chemical potential, leading 
to fj ~ {7r/fl)B{p). 

We now show that B{i') is actually proportional to the 
energy- resolved diffusivity of states at energy v. In the 
case of independent electrons which is of interest here, 
the electron mobility can be evaluated at any finite den- 
sity via Eq. (151, using the following expression for the 



(14) compressibility 



dn 
dp. 



dv 



dv 



p{y). 



(20) 



From Eqs. (15 1 and (191, and defining the diffusivity of 

(21) 



states at energy v as 

D{v)^{n/e)[B{v)/p{v)] 
we can rewrite the mobility as 



Me 



^dv {-l£)p{^) D{v) 



kuT 



Jd. 



(22) 



Our aim is to address the charge dynamics in systems 
where localized and itinerant states coexist in different 



which has explicitly the form of an average over 
energy with a probability distribution Wlv) = 



4 



^— p{v)/ J dv p{^)- In the limit of vanishing 

density, by taking the /i — > — cx) hniit appropriate to a 
non-degenerate semiconductor in Eq. (22 1, we find 



e Sdvp{v)e-P'^D{v) 



III. APPROXIMATION SCHEMES 

A. Semiclassical Kubo bubble approximation 

A powerful approximation scheme to calculate the 
carrier mobility is to evaluate the Kubo formula us- 
ing the exact electron propagators obtained in the limit 
of static molecular displacements, but neglecting vertex 
corrections^'^. Evaluating the single-particle propagators 
in the static limit is justified in virtue of the low fre- 
quencies of the intermolecular vibrations that couple to 
the electronic motion. On the other hand, the neglect 
of vertex corrections amounts to dropping the quantum 
interference proc esses that are responsible for Anderson 
localization,EnEll in the spirit of the semiclassical approx- 
imation. It corresponds to replacing the function Biy) 
appearing in Eq. ( 22 ) by the factorized expression 



B{v) = tv ( {p(y)) J {p{v)) J ) , (24) 



where (...) means an average over disorder variables 
(the averaging procedure will be defined in the following 
Section). In diagrammatic terms, only the elementary 
particle-hole "bubble" — a convolution of two spectral 
functions — is retained in the evaluation of the current- 
current correlation function. 



While Eq. (24) neglects particle-hole quantum corre- 



lations, it still accounts for non-trivial interaction effects 
contained in the single-particle propagators, which are 
calculated exactly. In particular, it is able to capture 
those aspects of the transport mechanism which stem 
from the dual nature of the electron states. It is there- 
fore superior to the usual Bloch-Boltzmann treatment in 
that it can account for both the coherent motion of band 
states and the incoherent motion of tail states. The con- 



volution Eq. ( 24 1 is actually analogous to the form that 



applies in the limit of large lattice connectivity under- 
lying dynamical mean field theory, and that has proven 
successful to address the crossover from band-like to hop- 
ping motion of small polarons non-perturbatively.'^^ It 
reduces to the static treatment of the gaussian disor- 
der model presented in Refs. [33] and [34| in the classical 
limit where the electron bandwidth is neglected, which is 
appropriate for narrow-band amorphous semiconductors 
and polymers in the strong disorder regime. 



B. Relaxation time approximation (RTA) 

The theoretical framework developed in Sec. II allows 
us to restore the backscattering processes leading to An- 
derson localization, i.e. those that are neglected in Eq. 



(24 1, by performing a phy sically t ransparent relaxation 
time approximation 

(RTA)P2l23l^ The idea underlying 
the RTA is to express the dynamical properties of the 
electronic system under study in terms of those of a suit- 
ably defined reference system from which it decays over 
time, and that can be solved at reasonable cost. In the 
semiclassical theory of electron transport, for example, 
one starts from a perfectly periodic crystal and describes 
via the RTA the decay of momentum states due to the 
scattering by impurities or phonons. The idea here is to 
find an alternative reference system to start with, so that 
quantum localization effects are built-in from the begin- 
ning. We now show that this can be achieved by starting 
from the exact description of a "parent" localized system 
where the disorder is assumed to be static. The RTA can 
then be used to restore the disorder dynamics related to 
the low- frequency lattice vibrations.^'' ^^l 

Let us consider an organic semiconductor where the 
disorder variables (i.e. the molecular positions) fiuctuate 
in time over a typical timescale Tin. At times t <^ Tim the 
molecular lattice appears to the moving electrons as an 
essentially frozen disordered landscape. The velocity cor- 
relation function C+{t) [cf. Eq. ([2|] then coincides with 
what would be obtained if the disorder were static, which 
we denote C^'^(t) (our reference system). In particular, 
the buildup of quantum interferences underlying Ander- 
son localization — which occurs on the scale of the elastic 
scattering time t^i — is realized provided that Tei < Tm- 
In this time range, the organic semiconductor therefore 
exhibits all the features of a truly localized electronic 
system. Quantum interferences that were present in the 
parent localized system are instead destroyed at longer 
times because, due to the lattice dynamics, the electrons 
encounter different disorder landscapes when moving in 
the forward and backward directionJ^. The form 



C+(t) = C^°'=(t)e-*/" 



(25) 



is the simplest form that is able to capture such decay 
process. Transforming Eq. (25) to Laplace space, results 
inC+(p) = C^°^(p+l/r„). 

The corresponding diffusion coefficient can be straight- 
forwardly obtained from Eq. (11). We see that starting 



from a localized system with a vanishing diffusion coeffi- 
cient, Ci|°'^(p — > 0) = 0, the RTA restores a finite diffusion 
coefficient which equals the diffusivity of the localized 
system at a time Tin- This result can be expressed as 



2t- 



(26) 



which is analogous to the Thouless diffusivity of Ander- 

-2 



son insulatora^^. Correspondingly, the quantity Lf^^{Tin) 
evaluated through Eq. ([7| for the reference localized 
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system acquires the meaning of a transient localization 
length for the actual dynamical system, as it represents 
the typical electron spread achieved after the initial local- 
ization stage, and before diffusion sets back in at t > Tin 
(see Figs. 1 and 3 in Ref. [13] as well as Fig. 1 be- 
low for a real-time illustration of this behavior). The 
emerging physical picture of the electronic motion that 
follows from the RTA Eq. (261 is quite different from 



the usual semiclassical picture, where disorder and lat- 
tice vibrations cause rare scattering events on extended 
electronic states. The present scenario rather describes 
electrons that are prone to localization but can take ad- 
vantage of the dynamics of disorder to diffuse freely over 
a distance Lioc{Tin), with a trial rate 1/Tj„. As will 
be shown in Sec. V, the RTA essentially reproduces 
the results obtained from more tim e-consuming mixed 
quantum-classical simulationJ^ilSIini^ and is free from the 
known drawbacks of these approaches. 

Before presenting model-specific results in Section [V] 
we analyze in more detail how the energy-resolved quan- 
tities of Sec. II C translate into the RTA language. From 
Eqs. (15) and (26) the RTA mobility in the low density 
limit is 



nksTQ 2Ti, 



(27) 



In the spirit of Eq. ( 23 1 , the transient localization length 
Lioc{Tin) can be expressed in terms of its energy resolved 
equivalent, £{t,i/), i.e. the spread reached by electronic 
states of energy v at time Ti„ , as 



li^^^ LUnn) ^ /d^p(^)e-/^-^^(r„„^) 



[see Appendix [c] for an explicit expression of (^{rin, v)^^- 
Combining Eqs. (p3|), (27) and (28) we recognize the 



energy-resolved diffusivity 



It- 



(29) 



which relates directly the conduction properties of the 
electronic states to their localization length in the parent 
localized system. 

Finally, from the considerations of the preceding Sec- 
tion we can derive the following relation: 



(30) 



Eq. (30) expresses the electron mobility in the RTA in 



terms of the optical conductivity of the reference local- 
ized system, whose mobility strictly vanishes. This result 
deserves a few comments. From scaling theories of local- 
ization, a finite d.c. conductivity is customarily obtained 
by taking the optical conductivity to saturate at a cut- 
off frequency of the order of the inverse of the inelastic 
scattering time, ct^.c. = c^ioc{uj ~ r"^). In Eq. (30 1, 



instead, the inelastic scattering time enters into the de- 
termination of the mobility via a weighted integral (i.e. 
through a lorentzian convolution) that involves the con- 
ductivity at all frequencies. The mobility Eq. (30 1 can 



therefore be quite different from the value obtained from 
the usual thumbrule. Our scheme is also conceptually 
different from the approach used in Ref.'^Sl. There the 
mobility was obtained by performing a Lorentzian con- 
volution of the optical conductivity itself, with a phe- 
nomenological broadening F = l/ri„ that was assumed 
to originate from the quantum fluctuations of the molecu- 
lar vibrations instead of the classical molecular motions. 
Apart from its different physical content, the method of 
Ref. [36] provides, for a given value of r"^^, a lower esti- 
mate for the mobility than Eq. ( 30 ) . 



IV. MODEL AND METHOD 

We now apply the theoretical framework developed 
in the preceding Sections to a model relevant to or- 
ganic semiconductors, that accounts for both the in- 
trinsic dynamica l disorder caused by inter-molecular 
motionjnnSminni a,nd the fluctuations of the molecular 
site energies that are assumed to originate from extrin- 
sic sources disorder. Specifically, we consider the follow- 
ing tight-binding Hamiltonian, for electrons or holes on 
a one-dimensional molecular lattice 



H 



(ij) 



cfco + h.c. 



(31) 



where are molecular site energies, and tij are inter- 
molecular transfer integrals between nearest neighboring 
molecules. In a perfect crystal all site energies are 
equal and can be set to zero without loss of generality. 
Static disorder leads to variations of the site energies with 
a statistical distribution P{ei). We are interested here in 
the effects of energetically distributed disorder, as op- 
posed to disorder centers with a definite energy (such as 
specific point defects). Correspondingly, we take P{ci) 
to be a gaussian of variance A as a representative case 
study. In addition, the coupling of the electrons to the 
vibrations of the molecules induces a dynamical disorder 
in the inter-molecular transfer integrals tij which depend 
on the molecular positions Ri,Rj. These fluctuate on 
a timescale governed by the relev ant vibr ational modes, 
whose frequency we denote as oj^^^^^. We assume a 
linear dependence of tij on the intermolecular distance, 
t,j = toil - a{R, - R,)]. 

One-particle properties — i.e. the properties that de- 
rive from the electron Green's function, such as the spec- 
tral function, the quasiparticle lifetime or the density of 
states (DOS) — can be efficiently evaluated by treat- 
ing the molecular degrees of freedom as static. This is 
justified because the frequencies of the inter-molecular 
vibrations that couple to the electron motion are much 
smaller than the band energy scale, which follows from 
the large molecular mass. For example in rubrend^ 
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Wo ^ 4-9meV, to =i ISOmeV, so that < (seetnOHini 
for different compounds). The static approach amounts 
to treating the positions Ri as classical variables dis- 
tributed according to a Gaussian distribution of ther- 
mal origin" " FWi?^) cx ex-p{-Muj^Rf /2kBT) (M is the 
molecular mass) !^ In practice, a numerical solution for 
the electronic problem is obtained for each given con- 
figuration of {Ri} and {ei} and then averaged over the 
disorder configurations. The electronic properties of the 
model Eq. ( [3l| ) in the static limit depend on two di- 
mensionless coupling parameters: A/tg, that controls the 
amount of extrinsic disorder, and A = a^to/ {2A'Iluq), the 
electron-molecular lattice coupling parameter. From the 
latter, the variance of the intrinsic thermal fluctuations 
of the inter-molecular transfer intergrals is obtained as 
s = VSAT^o. 

The static limit described above leads to a strictly 
vanishing particle diffusivity for all states. This can be 
seen from Eq. (29), which vanishes when Tin — ^ oo. 



The dynamical nature of the inter-molecular vibrations 
must therefore be accounted for in order to address the 
transport properties of electrons. Both approximation 
schemes described in Sec. |III| accomplish this task by 
taking the solution of the static model as a starting point. 
The technical details are described in what follows. 

(i) To obtain the spectral function p^v) needed in 
the semiclassical Kubo bubble (KB) approximation we 
adopt an algorithm based on regularization of the tri- 
diagonal recursion formulas for the electron propagator 
(see Ref.ti3]for details). Using this method, system sizes 
up to iV = 2^^ sites can be achieved. The spectral func- 
tion is then obtained after averaging up to 6- 10^ different 
realizations of the Ri and e^. The mobility is directly ob- 
tained using Eq. (24). The average diffusivity C+(p) is 



obtained by evaluating the optical conductivity via Eq. 
(A20) applying the same factorization as in Eq. (24 1, 
and then using Eqs. ([8| and (10). 

(ii) To evaluate the average diffusivity C^"'^{p) needed 
in the RTA we use standard exact diagonalization tech- 
niques on chains of up to = 2^ sites. The functions 
of interest are then calculated via their Lehman repre- 
sentation as a sum over the resulting eigenstates (see 
Appendix A). Since we are considering electrons mov- 
ing in the time-dependent potential of the fluctuating 
molecular lattice, it is a natural choice to associate the 
scale Tin in the RTA Eq. ( 25 ) with the typical timescale 
of inter-molecular vibrations. Assuming a single vibra- 
tional mode with a frequency cjq — bmeV , in the range of 
the relevant inter-molecular vibrations in rubrene, we set 
Tin = 1/wo = 10^^"^s. This assumption has been shown 
in Ref. to be consistent with t he op tical absorption 
data available in Rubrene OFET^^. This choice is 
also consistent with the results of dynamical Ehrenfest 
simulations in Ref. [UJ as it correctly reproduces the 
departure from localization observed at time r ~ I/wq. 
However, a more precise estimate of Tin from the Ehren- 
fest results is prevented due to the inaccuracy of this 
method in the long time limit (see below). 



(iii) For comparison, we shall also present results ob- 
tained with the method of Refs.^^ (termed Ehren- 
fest in the following). The dynamics of the molecular 
positions Ri are then included explicitly by adding a vi- 

2^ to the Hamilto- 



brational term Hyu, = 



nian Eq. (31 1 {Pi are the conjugate momenta of the Ri) 



The electron diffusion is obtained via Ehrenfest quantum- 
classical dynamical simulations: the Ri are treated as 
classical variables subject to forces evaluated as averages 
over the electronic state obtained fr om th e solution of its 
time-dependent Schodinger equatiorP^J^^. We average up 
to 12800 initial conditions on a 1024-site chain, with the 
initial molecular displacements and velocities taken from 
the corresponding thermal distribution. For each initial 
condition we use a different set of disorder variables e^. 

(iv) By artificially freezing the Ri variables in the sim- 
ulation we obtain a formulation in the time domain of the 
static problem described at point (ii), that we shall refer 
to as Ehrenfest-S. From Eqs. ([T]) and ([8| the calculation 
of the time dependent mean square displacement of Eq. 
^ can then be used to obtain the quantities Lf^^{T) and 
C+(p) from a Laplace transform. 



V. RESULTS IN THE LOW DENSITY LIMIT 

A. Time-resolved diffusivity and transient 
localization 

Fig. [l] illustrates the function C+(p) defined in Eqs. 
([6])-([8]), as calculated from the different methods out- 
lined at the end of the preceding Section. This quan- 
tity has the meaning of a diffusivity averaged up to a 
time T — 1/p and therefore provides direct information 
on the quantum dynamics of electrons as a function of 
time. We fix A = 0.17, which is representative for the in- 
trinsic electron-vibration coupling in rubren^i^l^ and set 
the temperature to T = 0.2to = 300K. The numerical 
results for a single electron in a pure organic semicon- 
ductor (A — 0) and in the presence of extrinsic disorder 
(A = O.bto — GbmeV) are reported in Figs. [T^ and [Tja 
respectively. This value of A is close to the one that was 
derived from the analysis of angle-resolved photoemission 
spectra (ARPES) in crystalline pentacene filmpSl. 

All the methods outlined above yield a ballistic time 
evolution in the short-time limit. It can be shown that 
the average diffusivity at short times obeys exactly C+ = 
2{v'^)/p — 2(u^)r, which is ruled by the average band 
velocity (w^) (cf. Appendix |b|. The ballistic regime is 
followed by a flattening of the average difusivity due to 
the onset of scattering processes, occurring on a timescale 
which we identify with the elastic scattering time, Tej^. 
From Fig. [l] we estimate approximately Tei ^ 10~^'*s 
in the pure case (A — 0), and t^i ^ 5 ■ 10~^^s in the 
disordered case (A = 0.5to). 

In the long time limit, the different methods yield qual- 
itatively different behaviors reflecting the fundamentally 
distinct treatments of the intermolecular dynamics. Let 
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FIG. 1: Average diffusivity C+(p) defined in Eq. ([8| as 
a function of the timescale r = at a temperature 

T = 0.2fo = 3007^. Panels (a) and (b) are respectively with- 
out A = 0.0 and with A = O.Sto extrinsic disorder. The 
bold short-dashed curve (light blue) is obtained by solving the 
time-dependent Schrodinger equation on a 1024-sites chain 
(Ehrenfest-S). The dotted curve (magenta) is the result of 
the dynamical Ehrenfest approach with vibrational frequency 
a;o = 0.05to- The long-dashed line (green) is the RTA approx- 



imation with Tin = <^n 



The full (red) line is the result of the 
Open cir- 



(lOl 



semiclassical Kubo bubble approximation Eq. 
cles are results from exact diagonalization of a 256-sites chain. 
The arrows mark the elastic and inelastic scattering times (see 
text). 



us focus on the intrinsic case first, Fig. [iji. In the par- 
ent system with static disorder, the electrons are local- 
ized (Ehrenfest-S and ED, respectively blue short-dashed 
curve and open circles) . The existence of a finite localiza- 
tion length L(r) = const as t — >■ cx) implies through Eq. 
([S]) that the average diffusivity bends down and tends to 
C+ oc 1/t at long times. Restoring the dynamical nature 
of the inter-molecular vibrations via the RTA causes a 
departure from the localized behavior on the scale of the 
inelastic scattering time ri„, so that a diffusive behavior 



(C+ = const) is reached at long times (green long-dashed 
curve). We note that within the RTA the diffusion coef- 
ficient at r — oo is necessarily lower than the maximum 
attainable value, which is obtained when Tin — Td- 

Finally, the Ehrenfest method (purple, dotted line) 
also captures the departure from localization occurring 
at T ~ Tin- However, this method yields a spurious su- 
perdiffusive behavior at long times^^, which is testified 
by a marked upturn of the diffusivity. This drawback 
leads to an overestimate of the mobility, whose value can 
vary strongly depending on the chosen simulation time. 
For this re ason^ the m obilities obtained from Ehrenfest 
simulationJliniMllolli] should be taken with some care. 

The results reported in Fig. [T^, that have been ob- 
tained using microscopic parameters appropriate for pure 
rubrene (the organic semiconductor with the highest mo- 
bility reported to date), indicate that even in ideal sam- 
ples without extrinsic disorder the elastic scattering time 
at room temperature is shorter than the typical inelastic 
scattering time, t^i < Tin- This situation results from 
the combination of the large mass of the molecular units, 
which leads to low vibrational frequencies and therefore 
to large values of Tin, together with the typically weak 
inter-molecular transfer rates and their strong sensitivity 
to inter-molccular motions, which lead to short values 
of Tel- It is therefore expected to be a general feature 
of all organic semiconductors, resulting in a transport 
mechanism that is fundamentally different from that of 
inorganic materials. Specifically, a "transient" localiza- 
tion regime emerges at intermediate times, t^i < t < Tin, 
where the electrons tend to localize (the diffusivity bends 
down) before a constant diffusivity sets back in at r > rj„ 
(green long-dashed curve in Fig. [T^). When such tran- 
sient localization is realized, the diffusion coefficient at 
long times depends on the history of the system at this 
intermediate stage, being inversely proportional to the 
inelastic time Tin, cf. Eq. (26). 



The existence of a transient localization phenomenon 
invalidates in principle semiclassical treatments of elec- 
tron transport, which are instead successful in inorganic 
materials. To illustrate this point we show in Fig. [T^ 
the average diffusivity C+ obtained from the semiclas- 
sical Kubo bubble approximation (red full line). Be- 
cause in this approximation the backscattering processes 
at the origin of localization are neglected, the system 
evolves continuously from a ballistic to a diffusive behav- 
ior. Semiclassical approaches are therefore inadequate 
to describe electron transport in organic semiconductors 
where t^i < Tin- On the other hand, the present RTA is 
able to recover the semiclassical picture in the opposite 
regime where Tin < t^i- As discussed at the beginning 
of this Section we have that Ctf.{p) = 2{v^)/p at short 
times. Using the RTA Eq. ^ and Eq. ([u]) yields a 
diffusion coefficient D = {v'^)Tin, which is formally anal- 
ogous to the Bloch-Boltzmann result. 

The situation in the presence of extrinsic static disor- 
der is not qualitatively modified with respect to the pure 
case, as can be seen from Fig. lb. In particular, the in- 
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elastic scattering time remains unchanged, because it is 
determined by the intrinsic timescale of the intermolecu- 
lar vibrations. However, an increased amount of disorder 
shifts the onset of locaUzation Tei to shorter times. This 
enlarges the time interval where the transient localization 
phenomenon is effective, with a consequent reduction of 
the diffusion coefficient. 

Finally, one can wonder if the transient localiza- 
tion scenario, which was demonstrated here in the one- 
dimensional case, is robust in more realistic descriptions 
of OSC. For materials with a sizable in-plane anisotropy, 
i.e. such that the inter-molecular transfer rates in the 
transverse direction are much smaller than in the longi- 
tudinal direction, t± ^ to, the one-dimensional picture 
is expected to remain valid up to times t ~ h/tj_. This is 
the situation that applies to Rubrene, according to recent 
photoemission data.l^S Because the elastic scattering time 
is very short, of the order oi Ti/to or less (see Fig. 1), we 
have h/t± 3> Tel and nothing prevents the transient lo- 
calization to occur in this case. Backscattering processes 
remain relevant also in isotropic two-dimensional mate- 
rials. Although the timescale for two-dimensional weak 
localization is known to be longer than Tei ^ the present 
scenario should remain qualitatively correct also in that 
case due to the strong intrinsic disorder present in OSC. 



B. Temperature dependence of the mobility 

Fig. 2a shows the temperature dependence of the mo- 
bility as obtained from the RTA for different amounts 
of extrinsic disorder, A = 0, 0.2to and 0.5tQ. The low- 
est accessible temperature is set by the limits of validity 
of our classical treatment for the molecular vibrations, 
namely T ^ = O.OSto- The intrinsic mobility of pure 
compounds (A — 0, full red line) is a monotonically de- 
creasing function of temperature, even though the mi- 
croscopic transport mechanism is far from a conventional 
band transport, as discussed in Sec. |V A[ Depending on 
the explored temperature window, a power law temper- 
ature dependence, T~'^, can be identified. In practice 
the exponent 7 depends on how the transient localiza- 
tion length in Eq. (27) varies with temperature, due to 



the thermal increase of inter-molecular disorder. 

The behavior of L'iod'^in) is shown in the inset of 
Fig. [2}i. In the temperature interval T = 200 - 600i^, 
the transient localization length decreases steadilj'^ as 
^fociZia) ^ a^^o/(AT). Substituting t his ex pression in 
Eq. ([27) leads to a mobility fie oc T"^ EE^. Moreover, 
using explicitly the definition of A given in Sec. IV and 
the relation l/ri„ ~ loq we obtain that the mobility in 
this regime increases with the third power of wq, while 
it is independent of the transfer integral Iq. An analo- 
gous calculation in the semiclassical regimd^ yields for 

the model under study a mobility proportional to uiq and 
1/2 

to • In both cases, increasing ujq results in an increase 
of the charge mobility, because it suppresses the effects 
of dynamical inter-molecular disorder. This observation 
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FIG. 2: Mobility as a function of temperature calculated (a) 
within the RTA and (b) within the Kubo bubble approxima- 
tion. The inset of panel (a) shows the behavior of the tran- 
sient localization length Lf^^{Tin) with l/ri„ — ujq = O.OSto. 
Mobilities are expressed in units of fio = ea^/h, with a the 
lattice spacing (/xq = Icrn} jVs taking a = 7.2A for rubrene). 
In this and all subsequent figures, we take to as the unit of 
energy. 



suggests a possible strategy for the design of high mobil- 
ity materials: rather than optimizing the inter-molecular 
TT overlaps that control the value of to, it could be advan- 
tageous to stiffen the inter-molecular vibrations either via 
an appropriate tayloring of the inter-molecular structure 
(e.g. by molecular functionalization) or via the interac- 
tion with a substrate (as in self-assembled monolayers). 

The localization length becomes a weakly increasing 
function of T at very high temperatures, where a vibra- 
tionally assisted electron motion arises via the large fluc- 
tuations of the inter-molecular distances. This results in 
a weaker power law dependence of the mobility with ex- 
ponent 7 < 1, which was termed "mobility saturation" 
in Rcf.SI. Although such high temperatures are not at- 
tainable experimentally in rubrene, the mobility satura- 
tion regime could actually be observed in materials with 
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a stronger electron-vibration coupling constant A than 
that considered here. 

The inclusion of extrinsic disorder (A 7^ 0) causes a 
downturn of the mobility at low temperatures, that is 
reminiscent of a thermally activated behavior, i.e. fie 
increases with T. This behavior reflects a crossover be- 
tween the extremely short Lf^^^T^n) obtained at low tem- 
peratures and the larger intrinsic value at higher tem- 
peratures, as shown in the inset of Fig. (2^.1^ The loca- 
tion of the crossover from extrinsic to intrinsic transport 
depends on the amount of extrinsic disorder, so that it 
can vary experimentally depending on the material and 
device quality. Correspondingly, a variety of behaviors 
ranging from thermal activation to a power law decrease 
can be realized in the experimental temperature window, 
which is possibly at the origin of the different tempera- 
ture dependent mobilities observed in organic FETs. 

The results of the semiclassical Kubo bubble approx- 
imation are shown in Fig. for comparison. Despite 
the profound differences in the two descriptions of charge 
transport, the overall behavior obtained in the explored 
temperature interval is qualitatively similar. From a fun- 
damental viewpoint, the fact that a thermally activated 
behavior is obtained for A 7^ in both the RTA and the 
Kubo bubble approach indicates that the corresponding 
hopping processes (incoherent jumps from molecule to 
molecule) are already present at the semiclassical level, 
being captured by the Kubo bubble approximation in the 
strong disorder limit .'^ This result means that the hop- 
ping behavior is not related to the quantum (backscat- 
tering) localization corrections, as the latter are not con- 
tained in the semiclassical treatment. As we proceed to 
show, the crossover from the intrinsic to the thermally 
activated regime can be explained in terms of the com- 
petition between highly conducting states located in the 
band range and weakly mobile states located in the band 
tails — a competition that is captured by both the RTA 
and the Kubo bubble approximation. 



Energy resolved diffusivity and localization 
length 



Based on Eq. (231, the electron mobility in a non- 



degenerate semiconductor arises from a weighted average 
of the energy-resolved diffusivity D(v) via the thermal 
population of electronic states. The latter is measured 
by the thermally weighted DOS, 



(32) 



which represents the normalized probability of occupa- 
tion of the states at energy ly. The functions p{iy), D{v) 
and W{v) are analyzed below. 

Fig. [3] illustrates the evolution of the DOS p{v) as a 
function of increasing thermal disorder in a pure sam- 
ple (a) and upon increasing extrinsic disorder (b). The 
DOS of a perfectly ordered crystal is shown for reference 



(dashed). In both cases, the Van- Hove singularity mark- 
ing the edge of the one-dimensional band at v = —2tQ is 
rounded off and shifts deeper in energy, indicating an in- 
crease of the effective bandwidth.'-'^^''*^' In addition, tails 
are generated beyond the range of band states. Both 
the bandwidth increase and the extension of band tails 
are controlled by the amount of disorder. This can be 
quantified through the variance s = VSATto of the ther- 
mal fluctuations of the inter-molecular transfer intergrals 
in the intrinsic cas^^^ (Figjs^), and by the spread A of 
molecular energy levels in the extrinsic cas^^ (FigjsjD) . 

The diffusivity D{i') obtained through the RTA in the 
pure case (A = 0) is shown in Fig. (right axis scale). 
The diffusivity is directly related, through Eq. (291, to 
the square of the energy-resolved transient localization 
length, i'^{Tin, v), which is also shown on the same figure 
on the left axis. Analogous estimates for the localization 
length in the static disorder limit — >■ 00 can be found 
in Refs. I13I21I The comparison with the DOS of Fig. [3}i 
allows us to identify two distinct regions in the electronic 
spectrum, separated by a crossover region of width ~ s 
around the band edge. States located in the band region 
have a large diffusivity, that is strongly suppressed upon 
increasing the thermal disorder. Tail states induced by 
disorder below the band edge instead have a much lower 
diffusivity as a consequence of their more localized char- 
acter. The diffusivity of tail states is essentially tempera- 
ture independent and corresponds in our one-dimensional 
model to a minimum localization length of approximately 
one lattice spacing, ^(Tj„,i/) « a. We note that the ex- 
istence of two distinct characteristic values of the local- 
ization length is in agreement with recent ESR measure- 
ments performed on pentacene transistors!'*^''^^ 

The relative importance of band and tail states in 
the transport mechanism is determined by the weight- 
ing function W{v), which is shown in Fig. |4|d. At tem- 
peratures T ^ 0.5io, which includes the experimentally 
accessible range, the function W{v) is peaked right in the 
crossover range that separates the band and tail states 
(cf. Fig. [3]), with a sizable overlap on both sides. As ex- 
pected for two conduction channels in parallel, the elec- 
tronic transport in this case is dominated by the channel 
whose diffusivity is largest, i.e. the band states. Corre- 
spondingly, the temperature dependence of the mobility 
is governed by the suppression of the diffusivity in the 
band range, that is illustrated in Fig. |4^. Upon increas- 
ing the temperature, the weighting function progressively 
broadens and shifts towards the tail states. These even- 
tually become the dominant transport channel, leading 
to the mobility saturation observed in Fig. [2^. 

The energy-resolved diffusivity and transient localiza- 
tion length obtained in the presence of extrinsic disorder 
(A = O.Sio) are illustrated in Fig. |4]:. The main dif- 
ference with the pure case shown in Fig. is that the 
crossover region separating tail and band states is now 
broadened by an amount cx A. Nevertheless, the typi- 
cal values of the diffusivity both in the band region and 
in the tails remain close to their intrinsic values, pro- 
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FIG. 3: (a) Density of states in a pure crystal with intrinsic inter-molecular disorder and A = 0. The different curves correspond 
to increasing temperatures from T — O.lto to T = 0.7to- The dashed line is the DOS of the perfect crystal in the absence of 
thermal disorder (A = 0). (b) DOS at T = 0.2to for several values of the extrinsic disorder A. 
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FIG. 4: (a) Energy-resolved transient localization length squared t^ijin^v) (left axis) and diffusivity D{y') (right axis) in a 
pure OSC, for increasing temperatures from T = O.lto to T = O.Qt p. P arameters are the same as in the preceding figures. The 
length unit is the lattice spacing a. (b) Weighted DOS [see Eq. |32l] at the same temperatures, (c) and (d) show the same 
quantities in the presence of extrinsic disorder (A = 0.5). In panel (d) the temperature T = 0.15to has been added to highlight 
the crossover from extrinsic to intrinsic transport. 



vided that the temperature is not too low. This is more 
clearly seen in Fig. [5^, which shows the RTA diffusiv- 
ity at increasing values of A for fixed T = 0.2io- At 
lower temperatures (T < O.l^o) the extrinsic disorder 
eventually becomes dominant and the results tend to re- 
cover those obtained in the absence of intrinsic electron- 
vibration coupling (A = 0, dashed curve in Fig. |4]:). 

Since for A = O.Sio the diffusivity is a monotonically 
decreasing function of T for all states, the origin of the ac- 
tivated behavior of the mobility observed at low temper- 



atures in Fig. [2^ has to be sought elsewhere, i.e. in the 
weighting function W{v). As illustrated in Fig. [4|i, the 
behavior of W(v) in the presence of extrinsic disorder is 
opposite to that of the pure case. Fig. |4]d: at low temper- 
ature (here T = Clio) the peak of the weighting function 
is located deep in the tail states and it moves towards the 
band upon increasing the temperature. Tail states with 
a low diffusivity therefore dominate the transport mech- 
anism at low temperature, while the intrinsic regime is 
progressively recovered upon increasing the temperature. 
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FIG. 5: Energy resolved localization length squared {Tin , i') 
within the RTA (a) and Kubo bubble approximations (b). In 
both cases we take l/rin = ujq — 0.05to 



The crossover between these two regimes is signaled by 
a maximum in the mobility of Fig. [2^ at a temperature 
that we denote T* . By comparing the variances of intrin- 
sic and extrinsic disorder given at the beginning of this 
Section, i.e. setting s ~ A, we obtain the following es- 
timate for the crossover temperature: T* ~ A^/(8Ato)- 
Taking A = 0.5^0 and A = 0.17 gives T* = O.lSto, in 
good agreement with the data of Fig. [2^. The predicted 
crossover temperature for A — 0.2io is T* = O.OS^o, out- 
side the studied range. 

The Kubo bubble approximation yields qualitatively 
similar results for the temperature and A dependence 
of the diffusivity 0(1/) (of course the DOS and weighting 
function W{i') are exactly the same, as they are obtained 
in the common static limit). The main differences be- 
tween the two methods are quantitative and arise from 
the inclusion or not of backscattering effects, i.e. they 
are indicative of the relevance of vertex corrections in 
the disordered system under study. As shown in Fig. [5] 
the diffusivity of band states in the Kubo bubble approx- 
imation is larger than in the RTA, corresponding to the 
fact that for band states the Kubo bubble essentially re- 
covers the Boltzmann transport theory^'^ where quantum 
localization phenomena are absent. This leads to larger 
values of the intrinsic mobility than in the RTA, as seen 
in Fig. [2]3. Concerning tail states, the opposite is true. 
In the RTA the localization length and hence the diffu- 
sivity appear to be bound from below, which provides a 
lower bound to the mobility at low temperatures: taking 



FIG. 6: Mobility as a function of carrier concentration per 
molecule within (a) the RTA and (b) Kubo bubble approx- 
imation, for T = 0.2to- The inset shows the variation of 
mobility with respect to its zero density limit. 



i{Tin,i^) > a from Fig. [5^ yields fie > ea? / {2TinkBT). 
This lower bound is absent in the Kubo bubble approx- 
imation, where D(y) vanishes asymptotically for nega- 
tive energies, being itself proportional to the DOS. It 
can actually be showrP^l that a behavior of the form 
obtained from the Kubo bubble in the 
limit of strong disorder, A ^ to. As a result the ther- 
mal activation at low temperatures is more pronounced 
in the Kubo bubble approximation than in the RTA. The 
two effects discussed here are at the origin of the much 
sharper crossover from extrinsic to intrinsic transport ob- 
tained in the Kubo bubble compared to the RTA, cf. Fig. 



m 



VI. 



DENSITY DEPENDENCE OF THE 
MOBILITY 



The mobility obtained through Eq. (22) at finite 
electron concentration is shown in Fig. [6j The pa- 
rameters are the same as in Fig. [l] i.e. A = 0.17, 
T = 0.2io = SOOiC. In both the RTA, Fig. |6}i, and 
Kubo bubble approximation. Fig. [6]d, we find a steady 
increase of the mobility with increasing density. This be- 
havior can be understood as arising from a progressive 
filling of tail states, which occurs through a shift of the 
chemical potential towards the band region as the elec- 
tron liquid becomes degenerate. This allows states with 
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FIG. 7: Mobility as a function of carrier concentration per 
molecule, obtained from the RTA for T = Q.2to and A = Q.2t 
(dashed line, same as Fig. 6a). The full line is obtained 
by replacing the definition Eq. ( |15[ ) by the expression ^.e = 
cr/(ne). 



a higher diffusivity to be populated, via a shift of the the 
factor —df/dv in Eq. (22 1. This argument is only quali- 



tatively correct however, because it neglects the fact that 
the diffusivity D{v) itself depends on the density, which is 
instead correctly incuded in the results of Fig. [6j Based 
on the same argument, increasing the density of carriers 
will shift the crossover between the extrinsic and intrin- 
sic regimes to lower temperatures. Such depinning effect 
can be achieved in OFETs through the application of a 
strong enough gate electric field. 

We see from Fig. |6]that the curves describing the den- 
sity dependence of the mobility for different degrees of 
extrinsic disorder are essentially parallel. The quantity 
jjLe — /ie(?T. = 0), whcrc the n = limiting value has 
been subtracted, is shown in the inset. Interestingly, 
it appears to be insensitive to the presence of extrin- 
sic sources of disorder: the curves show a clear collapse 
for different values of the static disorder, that persists up 
to fairly large carrier concentrations. This indicates that 
at the considered temperature, T — 0.2io — SOOiiT, the 
observed increase of mobility with density comes entirely 
from populating carriers with a strong band character 
(cf. Fig. [4]). The quantity /Xg — /ie(n — 0) could therefore 
be used to measure the mobility of band carriers even in 
samples with a sizable degree of disorder. 

Finally, because large carrier concentrations are now 
customarily obtained in OFETs via liquid gatin^SfilsiJ^ 
the correct definition of the mobility, Eq. (15), should 



be used to analyze experiments, instead of the usual ex- 
pression a jijie) that only holds in the limit of vanishing 
density. To illustrate this point, we compare the results 
obtained with the two definitions in Fig. [Tj Use of the 
low-density expression (black, full line) leads to an er- 
roneous result, as it incorrectly predicts a reduction of 
the mobility upon increasing the density. To understand 
this result, we observe that for degenerate carriers, the 
susceptibility Eq. ( 20 1 is given by the DOS at the chem- 
ical potential, i.e. on/d^i ~ p{p)- Using Eq. (15) we 
can then rewrite <j/{ne) ~ fiekBTp{p)/n, which is itself 



proportional to the DOS (cf. Fig. 3b). In the present 
model an increase of the density implies a reduction of 
the DOS at the chemical potential, explaining the behav- 
ior observed in Fig. [Tj 



VII. CONCLUSIONS 

Based on a theoretical formalism that relates the Kubo 
formula for the conductivity to the time-resolved diffu- 
sivity of electronic states, we have analyzed the electronic 
transport mechanism in a model that accounts for several 
key ingredients relevant to organic semiconductors: the 
existence of narrow electronic bands, the dynamical dis- 
order arising from the thermal vibrations of the molecules 
and the presence of extrinsic sources of disorder that are 
unavoidable in real samples and devices. 

The presence of strong dynamical disorder is intrin- 
sic to organic semiconductors and invalidates the usual 
semiclassical treatments of electronic transport that ap- 
ply to inorganic semiconductors, calling for a theoreti- 
cal approach that is able to treat quantum localization 
corrections in a controlled way. This is achieved here 
through a relaxation time approximation (RTA) that re- 
lates directly the carrier diffusivity to the localization 
properties of the electronic states. Within this theoreti- 
cal scenario, the deviations from semiclassical transport 
are understood as arising from a transient localization of 
electrons, that takes place before the onset of a true dif- 
fusive behavior at long times. This phenomenon appears 
to be a characteristic feature of organic semiconductors, 
where the typical timescale of inter-molecular vibrations 
is longer than the elastic scattering time. The transient 
localization scenario is supported by numerical simula- 
tions on the time-dependent diffusivitjff^ and by op tical 
conductivity measurements in Rubrene OFETtPSIill. 

Based on the present theory, the intrinsic transport 
mechanism in clean organic semiconductors is explained 
as the diffusive spread of localized wavefunctions rather 
than the scattering of delocalized waves by phonons and 
disorder. A power-law decay with temperature is pre- 
dicted for the intrinsic mobility, which results from the 
reduction of the transient localization length as the ther- 
mal disorder increases. Our results suggest that the in- 
trinsic mobility of organic semiconductors could be im- 
proved by tailoring crystal structures with stiffer inter- 
molecular bonds, as this would reduce the impact of ther- 
mal disorder on the charge transport. 

The inclusion of extrinsic disorder causes a crossover 
from the intrinsic power-law behavior, persisting at high 
temperature, towards a thermally activated behavior in- 
duced by carrier trapping at low temperature. Increas- 
ing the electron concentration induces a depinning from 
trapped states, leading to an increase of the mobility 
and a progressive suppression of the thermally activated 
regime. Our results for the concentration dependence 
of the mobility generalize the find ings obtained in the 
classical hopping limit to <C ApSEH to the high mobility 
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organic FETs of present interest, where the existence of 
electronic bands requires a quantum treatment of elec- 
tron motion. 

From a more general viewpoint, the present work 
demonstrates that the conductive properties of both pure 
and disordered organic semiconductors can be efficiently 
understood within a unified framework, by addressing the 
interplay between mobile states in the band region and 
strongly localized states in the band tails. The present 
results confirm and extend the considerations of Ref. [13] 
by allowing for a proper inclusion of quantum localization 
phenomena. Interestingly, the relationship between the 
energy-resolved properties of electronic states and the re- 
sulting mobility, that we have exploited here, could also 
be generalized to study how the intrinsic polarizability of 
the organic crystals affects the transport characteristics, 
as was recently proposed in Ref. [5] 



Appendix A: Detailed balance and spectral 
representation 



Expressing C+{t) in terms of its Fourier transform 
C+{t) = J ^e-*'^*C+(w) and using the reality of C+{t), 
which implies that i?eC+(w) and ImC^{uj) are respec- 
tively an even and an odd function of oj, we have 



C+ip) 



dijj pReC ^{u) + tj/mC_|_(a;) 



(A3) 



This equation can also be obtained using the analiticity 
of C+(w) in the complex upper half-plane with the help of 
Cauchy's residue theorem for complex integration. The 
two terms in Eq. (A3 1, respectively proportional to the 



real and imaginary part of C"*", bring the same contribu- 
tion to the integral as can be checked via the Lehman 
representation of the correlation function 



In this Appendix we set e = 1 in addition to ?i = 
fcs = 1. By introducing the Laplace transform of the the 
retarded current-current correlation function C+(i) 



C+{t) 



i^e-^^"| <n|J|m> |22cos(L^„,,„t), (A4) 



C+{p) = / dte-P'C+it), 



(Al) 



we can write the mean square displacement defined in 
Eq. as 



C+{P) 



where \n > and En are respectively the eigenvectors and 
eigenvalues of the Hamiltonian, that are supposed to be 
known, and aj„ ,„ = En — Em- The Fourier transform of 



(A2) Eq. (A4) reads 



ReC+{uj) = ^^e~^^"| < n|J|m > p[(5(w- w„,„j) +(5(a; + ( 



ImC+iio) ^ ^^e-^^"| <n|J|m> |2- 



2uj 



(A5) 
(A6) 



Using the result 
P 



1 



p 



(A7) 



we arrive at 



Jo 



^"^^i?eC+(c.). (A8) 



From this equation we can pro ve the limit used to obtain 
the Einstein's relation Eq. ([is]) , observing that asp ^ 0'^ 
^^2q^ <5(w) then limp_>o C'+ (p) = lim^^o C+{uj)/2. 

Eqs. ([3]) and Q are two relations which are not inde- 
pendent since the functions C+ and C_ are related by the 



detailed balance condition. The Fourier transforms of the 
retarded correlation functions C± (t) can be expressed as 
C±{uj) = C>(a;) ± C<(ti;) where C>(a;),C<(w) are the 
Fourier transforms 



C>(w) 



+ CXD 



dte'^-^+'^'>' < J{t)JiO) > (A9) 



+ 00 



C<(w) = / dte*("+"^)* < J(0)J(i) > (AlO) 
Jo 



and S an infinitesimal positive quantity. The detailed bal- 
ance condition reads^^ i?eC<(aj) = e~^"i?eC> (w), which 
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leads to 



ReC-{uj) = tanh ( ^ ) ReC+{uj). (All) 



which is Eq. ( 10 ) of the paper 



Applying the detailed balance Eq. (All) to Eq. (A3) 
yields 



C+ip) = 



duj 



2pReC-{uj) 



TT (p2 +a;2)tanh(;9w/2) 



(A12) 



and using the Kubo formula Eq. ^ with the definition 
Eq. (A2) we arrive at 



Jo 



(p2 +a;2)tanh(/3w/2) 



(A13) 



Similarly to Eq. (A5), the Fourier transform of the 



anticommutator correlation function C_ (w) can be ex- 
pressed through its Lehman representation as 



i?eC_(w) = ^^e ^^^"1 < n\J\m > ^ [(5(w + a;„,™) - - w„,™)] 



(A14) 



The two terms in the sum can be rearranged as 



< n\J\m > p(5(a; + - Em) 



(A15) 



We now consider single-particle Hamiltonians for which the number eigenstates \{na} > are such that H\{na} >= 
'^/3£0)|{'^q} >• Once expressed using this basis Eq. (A15) takes the form (in the grand canonical ensemble) 



ReC-{uj) = 7r(l - e 



-/3uj 



[{n}{rn} 



< {na}\J\{ma} > fS[uj + e„(7i„ - m„)] I . (A16) 



Since J is a single-particle operator, the individual eigenvalues obey Ua — rUa = Pa with Pa = —1, 0, 1. Thus in the 
sum appearing in the 6 function reduces to — e/3 where we have p-y = 1 when 7 = a, = — 1 when ^ — f3 and 
Py — elsewhere. The matrix element of J reads in this case 



a,l3 



(A17) 



where \a > are single-particle states. The grand canonical averages appearing in Eq. (AI6) can be performed leading 
to 

i?eC_(w) = 7r(l - e-^'^) ^ | < a| J|/3 > p < ?i„ > (1- < np >)5{uj + e„ - ep). (A18) 

Q,/3 



where we have made use of the vanishing of the diagonal elements of J. We can introduce a dummy integration 
variable by writing 



i?eC_(cj) =7r(l-e-^'^) 



/OO 
dvf{v)[\ - f{uj + i^)] ^ I < a\J\p > \^5{cj + v~ €p)5{v - e„). (A19) 
-00 „. o 



where /(i^) is the Fermi function. Identifying the diagonal part of the spectral operator 5{iy — e^) = —^Im < 
a|(i/ — H)^^\a > and taking into account that (1 — e~'^^)f{v){l — f{uj + j/)) = /(i^) — f{uj + v) we finally arrive at 



ReC_{oj) = 7r(l - e 



(A20) 
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Using Eq. (|A20|) and Eq. (|A12|) we obtain 



dv I duj 



1 [f{v)-f{uJ + y) 
a;2+p2 tanh(/3w/2) 



tr[p{v)Jp{v + cj) J] 



Taking the low density limit, 13^ — — oo, we obtain 



2 



Jo 



up- + jp- 



tr\p(y)Jp(y + J] 



(A21) 



(A22) 



where the normalization factor is Z = J diye ^'^p{v). 



Appendix B: Optical conductivity sum-rules 



Eq. (10) allows us to derive exact relationships be- 



tween the asymptotic expansion of C+ {p) [or equivalently 
L'^{p)] to certain integrals of the optical conductivity. Us- 
ing the definition Eq. (A2| we write Eq. (A13) as 



CM ^ r ^ 



du) 2puj 



a{uj) 



P+p^ tanh(/3a;/2)' 
Taking the formal expansion in powers of 1/p we get 

(-1)" 



(Bl) 



2n+l 



(B2) 



where 



Sn 



dcj 2^2" + ! 

TT tanh(/3w/2) 



Taking into account the definition of the Laplace trans- 
form we obtain 



C+ip) 



27i + ld"C+(i) 



n=0 



dt"-- 



(B4) 



The leading term of the asymptotic expansion [n = 0) 
gives C+{p) ~ C+it = 0)/p where C+{t = 0) =< .P >. 
Equating the coefficients of the expansions we get 



.0 = {-irs„ 



which reads for n = 

C+{t = 0) = 



did 



2uj 



TT tanh(^w/2) 



cr(w). 



(B5) 



(B6) 



Setting p = 1/r as in Fig. 1 we have that in the short 
time fimit C+(t) = C+{t = 0)t =< > r. While the 
RTA obeys this sum rule because the correlation function 
is exact at short times, the Kubo bubble approxi- 
mation does not, which results in the slight discrepancy 



observed in Fig. 1 in the ballistic regime. This can be 
related via Eq. (B6) to the different behaviour obtained 
for cr(w) in the two approximations and points to the 
relevance of vertex corrections at all frequencies. 

Appendix C: Boltzmann theory from the RTA 

The quantity i{p, v)"^ is evaluated in practice from 



e(p, v) 



did 



1 



p{v) Jq p2 ^_ ^2 tanh(/3w/2) 
or, equivalently. 



tr[p{v)Jp{v+ijj)J] 
(CI) 



e{p,v) 



duj- 



1 



P{v) J-oo P^ + ^' 



;tr[p{v)j p{v + oj)j]. 



(C2) 

(B3) These relations are obtained from the definition Eq. 



(A22) via Eq. (28 I 



The Bloch-Boltzmann theory is customarily obtained 
from the RTA by taking the noninteracting band eigen- 
states as the reference system. Within the present for- 
malism this amounts to evaluating the trace in the above 
equation assuming that the eigenstates of the Hamilto- 
nian are good momentum eigenstates, |fc), with energy 
Ek- In that case J|fc) = Vk = dE^/dk and we obtain 



Using Eqs. (23) and (29), with t = 1/p the relaxation 
time for momentum eigenstates yields the Boltzmann 
form of the mobility: 



(rvl) 



(C4) 



where the thermal average is defined as {tvD 
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